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ABSTRACT. Survival models used in biomedical and reliability contexts 
typically involve data censoring, emd may also involve constraints in 
the form of ordered p>arameters. In addition, inferential interest 
often focuses on non-linear functions of natural model parameters. 
From a Bayesian statistical analysis perspective, these features 

combine to create difficult computational problems by seeming to 
require (multi-dimensional) numerical integrals over awkwardly defined 
regions. This pap>er illustrates how these apparent difficulties can be 
overcome, in both parametric' ^d non- parametric settings, by the Gibbs 
sampler approach to Bayesian computation. 



1 . Introduction 

For the Bayesian statistical analysis of other thEUi simple stylized 
models, the key tool for calculation is (multi-dimensional) numerical 
integration; see, for example. Smith et al (1987) for a review of 
available techniques. However, it is widely recognized that 

considerable numerical sophistication is typically required in applying 
these techniques, euid that this has thus far haunpered the development 
of routinely available, user-friendly, Bayesian computational methods. 

This is p>articularly true in the case of survival models used in 

biomedical and reliability contexts. Here, features such as data 
censoring, ordered parameters, eissumed convexity or concavity of dis- 
tributions, all conspire to produce complicatedly constrained regions 
over which numerical integrations eure required. Not surprisingly, the 
literature therefore contains very few instances of fully Bayesian 
analyses in survival contexts (i.e. , presenting full and accurate 
posterior summaries, rather than, say modal point estimates or second 
derivative uncertainty measures). 

Recently, however, Gelfand et al (1991) have shown that the Gibbs 
sampler approach to Bayesian computation (see, for example, Gelfand emd 
Smith, 1990, and Gelfand et al, 1990) effectively side-steps the 



seeming problems of awkwardly defined integration regions in truncated 
data and constrained parameter problems, and provides an eaisily 
implemented computational procedure. 

Our purpx)se in this paper is to illxistrate the simplicity and scope of 
the Gibbs saunpler for the routine Bayesian analysis of survival data, 
in both j>arajnetric and non-peu'ametric settings. In Section 2, we 
briefly review the Gibbs sampler and its general structure for 
constrained peu'ameter and censored data problems. In Section 3, we 
provide a remge of illustrations of how the methodology proceeds for a 
variety of peu'ametric models used in various survival modelling 

contexts. In Section 4, we give a non-parametric illustration of the 
methodology. 



2. The Gibbs Sampler, Constraints and Censoring 

2. 1 THE GIBBS SAMPLER 

In what follows, densities will be denoted generically by square 
brackets, so that Joint, conditional and m^Lrgi^al forms for 

random variables U, V appear, respectively, as It/, V], [t/|l^] and [1^]. 
Marginalization by integration is denoted by [I/] = /[l/jV] • [1^] . Given a 
collection of random variable^ with Joint density [t/^, 

shall refer to , s = 1, 2 k, as the full conditional 

densities. 

The Gibbs sampler is a simply described iterative stocheistic simulation 
scheme, whereby samples drawn from the full conditional densities are 
used to provide summaries of aspects of the Joint density. Given an 

arbiti'ary set of starting values, "?• “ random variate l/J is 

drawn from then a variate Is drawn from 

and so on until is drawn from tl/^|t/j ^k-1^' 

completes one iteration of the sampler and results In a generated 
1 1 

vector (U. t/.). Rep>eating this process, after i iterations we 

IK ^ . 

arrive at a generated vector (1/^ l/^). It can be shown that, under 

mild regulau'ity conditions (see, for example, Geman and Geman, 1984), 

as i-»« this random vector tends in distribution to a random vector 

having the Joint distribution It/j l/^] . 

One possible procedure for obtaining summau'ies of aspects of 

U/j l/j^l of interest is therefore the following. Run M independent 

p>arallel replications of the above seimpling procedure, so that, for 1 



the resulting generated vectors 



Judged to be sufficiently large. 



<‘'u "kj’ 

size M from 



J = 1.2 H 

'“i V- 



» can be regarded as an iid sample of 
Standaird moment, quantile or density 



estimation techniques can then be employed to estimate summairy features 
of interest. 



In the Bayesian inference context, is the Joint posterior 

density of unknown model parameters Univariate mairginal 

inference summaries for l/^, say, are simply obtained from 



. . . , . Generally, if marginal inference summaries are required 

for a specified function of (t/^, . . . ,f/^) , an iid sample of size M from 



the corresponding marginal density is immediately available on substi- 



tuting the generated vectors 




••V’ 



J ' 



functional form. Exploration of bivariate (or 
marginal summao'ies proceeds in an obvious manner, 
emd comments on the pragmatics of choices of i auid 
Smith (1990, 1991), Gelfand et al (1990, 1991). 



1,2 W, into the 

higher dimensional) 
For further detail, 
M, see Gelfand emd 



The Gibbs saunpler thus provides a simulation-baised alternative to 
direct numerical integration methods, and one which depends only on our 
capacity to generate random vairiates (reasonably efficiently) from the 
full conditional densities, We shall now look at this 

latter issue in the context of constrained parameter and censored data 
problems. Our discussion here will be kept to the minimum necessau^ to 
give the reader an appreciation of how the Gibbs sampler achieves 
crucial simplification. For a much more complete discussion, see 
Gelfand et al, (1991). 



2.2 MODELS WITH CONSTRAINED PARAMETERS 



Suppose a p«u'ametric model for data Y involves a k-dlmensional 

parameter vector 0, constrained to lie in a subset of R^. For 

simplicity of exposition, we shall aissume here that does not depend 
on y (as would be the case, for example, if some components of 0 were 
truncation parameters: see Gelfand et al, 1991). Suppose further that 
[y| 0 ], [ 0 l denote the (unconstrained) forms of likelihood emd prior, so 
that the (constrained) Joint posterior for 0 » (0j,...,0^) is given by 



[y|el tel ^ 
J [y|0l-(0l 



[0|y] 



Proceeding by direct numerical integration, we see that there is an 
Immediate problem in calculating the normalizing constant and subse- 
quent problems in performing (k-1 ) -dimensional integrals (over subsets 

of S^) to obtain meu'ginal density forms. 



However, consider now the full conditional forms required for the Gibbs 
seunpler. If S^(0j. J*i) denotes the cross-section of corresponding 
to the constraints on 6^ Imposed by for specified values of 6j, J*i, 
we have 



lejy, 6j, j*i] cc [y|e]-(e] 



e^es^(6j, j*i) 



Moreover, the constraints region for 6^ will typically be an interval, 
or a union of Intervals. 



It follows that the typical rauidom variate generation task required for 
the Gibbs sampler in this case, will simply be that of generating from 
specified univariate density shapes truncated to Intervals. This Is a 
relatively straightforward task: in any ceise, strikingly easier than 

high-dimensional numerical integration over complicated constraint 
volumes. ' , 



2.3 CENSORED DATA PROBLEMS 



Suppose a parametric model for data y = (y, y ) involves a 

1 n 

k-dimensional peurameter vector e, with likelihood defined by 



[y|0l = n (yje] . 



However, suppose that for J t m > 1 there exist Vj, 
instead of observing Yj exactly, we simply observe that 
so that the likelihood is actually given by 




such that, 

€ iVj. Wj). 



m-1 ^ r J 

n [y |e] • n (y,|e] . 

i=l J=m ^Vj 

(We are here eissuming a simple, fully sj>ecifled censoring process for 
convenience of exposition. For a more general discussion, see Gelfand 
et aJ, 1991). 

In this case, a moment’s reflection reveals that the full conditional 
forms implied by the above likelihood combined with a prior (6] are 
not, in general, easy forms to sample from. In particular', the 



integral terms may have no closed-form analytic expressions, so that 
standard envelope rejection or ratio-of-unlforms seimpling techniques 
are not readily applicable. 

However, suppose we consider 7' = additional unknowns, 

so that the unknown model parameters are (e, 7'), with the data given 

by (7*, V, V). where 7" = (7, ^ V = and V = 

1 D-1 m n 

(.W ,...,1/ ). Consider now the full conditionals required for the Gibbs 
ID n 

sampler: 



[ej7*, V, W, 6j, J*i. 7') , 1=1 k , 

[7^ 1 7*, V, W, 6, 7^, s*r, s fc m] , r - m, ... ,n 

Careful examination of the conditioning variables reveals that the full 
conditionals for 0j 0^ reduce to 

[0^.|7, Oy J*i] , 1=1 k , 

the forms that would have obtained in the uncensored case! Typically, 
these forms present no difficulty for random variate generation. 

For the full conditionals for 7 , . . . , 7 , examination of the forms 

m n 

reveals that these reduce to 

W 

r ^ 

[7^|0] / J [7^|0] , r=m,...,n ; 

V 

r 

namely, the sampling distributions for the 7^ restricted to the ranges 
(Vj, Wj). Again, these typically present no difficulty for random 
variate generation. 

The trick of treating censored observations as unknowns in combination 
with the Gibbs sampler leads to simple Bayesian calculation strategies 
in otherwise intractable problems (see, also, Tanner and Wong, 1987, 
for a related manifestation of the idea). In the next section, we 
illustrate this concretely by detailing the forms of the Gibbs seunpler 
arising in a range of parametric models used in various kinds of 
survival studies. 



3. 



Illustrations For Parametric Survival Hodels 



3. 1 ORDERED BINOMIAL PARAMETERS 



Consider conditionally independent observations Binomial (n^, 

0^), i = 1,2 k , where it is known that 0^ s s...s 0^ and we 

seek to make inferences about the 0^ (or functions, thereof, such as 

0^^j-0^ or Problems of this kind arise, for example, 

in reliability development testing (Smith, 1977; Fard and Dietrich, 

1987), vrtiere stages 1 k correspond to successive improvements in 

reliability. 

If the Joint prior density is taken to be proportional to 



n 



^i-1 ^i-1 

(1-0J 



over the simplex, {(6j 0^ : 0 s 0^ s s. ..s 0^ 5 1}, by 

conjugacy the joint posterior [0|y] hzis the same form, with support S^, 
but with a^, replaced by + y^, - y^, respectively. 



Implementation of the Gibbs seunpler is now seen to be very simple. The 
full conditionals are given by 

[0^|y, 9j, J*iJ - Beta (a^+ y^, ~ > i * l,.. ,k , 

restricted to the interval 0. , s 0. a 0. . (0^ = 0, 0, , = 1), and 

1-1 1 1+1 0 k +1 

random vau'iate generation is straightforward. 

3.2 CENSORED REGRESSION DATA 



Schmee and Hahn (1979) modelled log-failure times of motorettes tested 
at four different temperatures by a straight-line regression of 
log-failure versus transformed temperatxire. Censoring occurred 

whenever a motorette had not failed at the end of the test period. The 

uncensored case likelihood (y|0] is assumed to derive from y. = a + 

2 

B^j + where ~ )V(0, <r), i=l,...,k, J = 1 n,, but the 



■ij 



actual data, Z, given by 



• - 



y 



if 



V. 

1 



Y s V 
ij i 



Y > V 
iJ i ’ 



where Is the total test time at temperature corresponding to X^. 

To Implement the Gibbs sampler, as Indicated In Section 2.3, we Include 
the unobserved Y.. (i.e. , those where Y.. > W ) as further iinknowns In 

2 

the model. In addition to the basic pzu'ameters of Interest, a,fi and <r . 

Given conjugate normal prior forms for a,p and an Inverse-gamma prior 
2 

for <r , It Is easily verified that the full conditional forms for a, ^ 
2 

and <r Eire straightforwardly Identified conjugate forms (normal, normal 
and Inverse-gamma, respectively) obtained as If all the Y were 

precisely observed. The full conditionals for the unobserved Y . . au'e 

2 

simply N(a tr ), restricted to the range > W^. Again, random 

variate generation from all these full conditionals Is unproblematic. 

3.3 TRUNCATED BIVARIATE NORMAL DATA 

Consider a blvEU'late normal process (X^, Y^), i= l,...,n, where some of 
the y. are not observed. One context In which such data arises Is In 

i 

paired survival time studies (using observed logarithms of svirvlval 
times), where observation (YJ of the second of the paired patients Is 

terminated when the first of the pair dies, so that Y^ Is observed only 

If y . s X . 

i i 

More precisely, we assume lid pairs (X^, Y^) such that for i = 1 n. 



X. 1 
1 


~ N - 




■®r 


• 


2 


®'l2 

p 




y. 

1 ■> 






•®2 . 




- ®'l2 


^2 J 





We observe the pairs (X^, Z^) with = Y^ If Y^ s X^; otherwise we 

observe (X^, •), where Z^ = • Indicates that Y^ > X^. Suppose that the 

prior for (0^, B^) Is taken to be bivariate normal with mean (p^, 

and covariance matrix V, eind that the prior for the covariance matrix, 

Z, say Is taken to be an Inverse-Wlshart, so that [Z = W <(pR) 
p) , with all the hyperparameters p^, Pg, V, p and R known. 

Interest focuses on meirglnal Inferences for 0^, 6^ and Z, but, 

following Section 2.3, unobserved values of Y^ are also treated as 

unknowns In specifying the Gibbs sampler. Defining = (X^, Y^), T = 



(T. r ). T = n"^r, +..+ r ). e = (e,. e„) and m » Cm,, m,). it is 

In In i<£ 

easily verified that 

[0|r, Z, Z] - NiVin~^Z + V)~^ T + n‘^Z(n"^Z + F)m. Cn^:"^ ♦ , 

iz"^|r, z. 0] = v{(zcr - 0)cr. - 0)' + pR)~^, n + p> 

i 



and 



D'JXj, Z.. 9. Z1 . W92 * <ri20-;2(Xj - 9,). . 

truncated to (X^, ») if Z^ = *, with degenerate at Z^ otherwise. 
The required random vau'iate generation is routine. 



3.4 WEIBULL PROPORTIONAL HAZARDS WITH CENSORING 

Consider a sxirvival time model in which the hazard function X(t; Z) , 
for an individual with covau'i^te values Z at time t, is given by 

X(t; X) = pt^~^ exp(ZR) , 



where P = O^, 0^,...,^^)' is a vector of unknown regression parameters 
and p > 0 is the unknown Weibull shape peu'ameter. 



If are explicitly observed survival times and t t are 

I n n+r m 



censored (T > t) lifetimes, 

Jth case, the likelihood is given by 



with Zj denoting covariate values for the 



TT 




7ft 

n exp[ -t^ e j 






. 



Cleau'ly, whatever the prior specification, the resulting 
(p 2) -dimensional posterior is awkward to handle using standerd 
numerical integration procedures. 



However, it is eaisily verified that the second partial derivatives of 
the log-likelihood with respect to each of the p + 2 unknown peu^eters 
are all non-positive (see Dellaportas and Smith, 1991). If the prior 

density is chosen to be log-concave, it follows that all the pxssterior 

full conditionals are log-concave. The import of this observation is 
that highly efficient methods exist for r^Lndom variate generation from 
log-concave densities (see, in particular, Gilks and Wild, 1991), so 



that routine, straightforward Bayesian calculation for widely used 
cases of proportional haizeu'ds models is possible (see Dellaportas and 
Smith, 1991, for wider exploitation of log-concavity). 



4. A Nonparametric Illustration 

4. 1 INTRODUCTION 

Nonpai'ametric Bayesian inference for the survival function with right 
censored data hzs been studied by Susarla and Van Ryzin (1976), and 
Ferguson and Phadia (1979). However, we often encounter the situation 
where some observations are censored from the left and some 
observations are censored from the right (see Turnbull, 1974, for 
references to papers addr'^ssing doubly censored data sets from a 
frequentist perspective). 

In this section, we study a nonparametric Bayesian approach to such 
problems, which allows us to incorporate prior beliefs and frees us 
from malcing a restrictive (p>€UTunetric) model assumption for the 
survival function. Specifically, we assume that the distribution 

function F of survival times has a prior given by Ferguson’s (1973) 
Dirichlet process, D[a) . The measure a can be written as NF^, where F^ 

is the prior meem of F and F^d ~Fq) / {N * 1) is the prior veriance of 

F. The larger N, the more strongly the prior specifies that F concen- 
trates around F^. 

In the doubly censored data case, it is very difficult to obtain an 

explicit expression for non-parametrlc Bayesian estimators even in the 
form of the pxjsterior mean. We shall show, however, that the Gibbs 
sampler approach, which augments the data by using latent variables 
that decompose the number of the censored observation into the possible 
numbers of observations falling into each interval, provides a 

straightforwardly computed numerical solution. As Illustrated in 
Section 2.3, this augmentation facilitates the sp>ecif ication of 

appropriate full conditional densities, particularly here for the 
survival functions given the latent variables. The iterated sampling 
scheme then allows us to approximate the posterior distribution of the 
survival function. 

4.2 THE MODEL 

He shall illustrate the approach lising a model similar to that studied 
by T\irnbull (1974), who proposed a self-consistent algorithm for 
computing the generalized maximum likelihood estimators. Here, we add 
the Dirichlet process prior to the model. 

Let r, , r„,...,T denote the true survival times of n individuals that 
12 n 



could be observed precisely if no censoring were present. The are 

independent euid identically distributed with distribution F; that is, 
F(.t) - P(T s t) for t t 0. We consider the ceise that not all are 

observed precisely. For each i, we assume that there are "windows" of 
observations and s V^) that are either fixed constants or 

random vau'iables independent of the {T^>. We observe 
= max I min(T^, wp, V^] . 

Moreover, for each item, we also know whether it is left-censored with 
Xj = V^, or right-censored with or a precisely observed time 

with X. = T.. 

i i 

We assume that items (or patients) are examined at discrete times (for 
exeunple, monthly) and that there is a natural discrete time scale 0 < 
t. < t_ <..., t , with observed deaths classified into one of the m 

id IB 

intervals (0, t,],(t,, t^]...,(t ,, t 1. Let 6. denote the number of 

lid m- 1 ."1 1 

precise observations (=) in the period denote the number 

of left-censored (s) entries at age t., and X. denote the number of 
right-censored (>) entries at t^. It is eissumed that the left-censored 
entries all occur at the end of eige period (t^, The data can 

then be summarized by the following tabulation: 



Type of obs. \ age (0, t^] (t^, 



(t 



m-1’ 






( = ) 
(:£) 
(>) 



m 



m 



m 



Let Pj * « 1 - denote the survival function 

tj, so that the likelihood function is prop>ortional to 



evaluated 



at 



j=i 

Let ^4 = P, , - P, for J = l,...,m and let 6 , = P . The prior 

J j-1 J j , , r 

process specifies that the distribution of the 0’s is the Dirichlet 
distribution 



where 



“j = "‘fo'V - • 



for J= 1 m * 1, with 



C = 



riN) 



rc; "■/ 



The posterior distribution of 0 = (0,, 0_, ...,0 ,0 , ) is known to be 

12 m flH-l 

a mixture of Dirichlet distributions (see Antoniak, 1974). In the next 
section, we show how the Gibbs saimpler side-steps the need for direct 
computation of this mixture. 



4.3 APPROXIMATION VIA THE GIBBS SAMPLER 



To employ the Gibbs sampler, we use the idea of Section 2.3 and 
introduce latent veu'iables that decompose the numbers of censored 
entries into the numbers of observations belonging to individual 
intervals. Let Z^j, . . . ,Zjj denote the random variables that count 

the number of observations in tij that might fall in the intervals (0, 

*2^ ^^j-r respectively, so that Hj = Z^j. 

Further, let Z..,,, ...,Z , . denote the number of observations in X. 

j+lJ m+lj j 

that might fall in the intervals ( t ,, t, (t ,, t ], (t ,oo], 

J 1 m-1 m m 

respectively, so that X . = 2,r 



Our objective is to summarize, via samples generated form the Gibbs 
sampler, the posterior distribution of 0 given the data. The posterior 
full conditional for 0 given the Z’s and the data, is easily seen to be 
an up-dated Dirichlet distribution depending only on the Z’s. The 
posterior full conditional for the Z’s given 0 and the data, is easily 
seen to be a product of multinomial distributions. Thus, suppose at 
the ith iteration step of the Gibbs sampler, we have the realization 

0-^ = (ej, ei ©Li), with 0^-1. We then up-date the Z 

1 2 DI+1 ", . J 



vairiables from the multinomial distributions as follows. For each 

J, J = l,...,m, we sample Z^*j^ , . . . ,Zj*j^ from the multinomial distribu- 



tion with sample size fij and peo'^eters rjj rjj, where = ej 

/ 1 J- Similarly, we seonple zj+J y» — ’^^1/’ 



J=1 



the multinomial distribution with sample size \j and p>arameters 

^j+lj ''m+lj’ ^jj “ ®j / ej for I = J + 1, ... ,m + 1. 

Having saimpled the Z random V 2 ur*iables, we then generate new 0 

variables from the Dirichlet distribution as follows. We compute, for 
each J, 




“j * *j 



02 



'2J 



and then saunple . . . , 6^tl) from the Dirichlet distribution 

1 02 02+1 

with parameters (yf*^ 

1 02+1 



By running M parallel independent replications of the sajnpler, aifter 
the ith iteration, we have ej , . and yf •...»yi^- » for 

Is, 2s 02+ 1 . S Is 02+ 1 . S 

s = The posterior distribution of 0^ for J = + 1 can 

then be approximated (for sufficiently lairge 1) by 



- H m+1 . 

- «- [ Beu(r I r^) , 

s=l k*l 

where Bet&{a.,b) denotes the beta density with parameters a and b. A 
posterior estimate of the 0^ is then given by 



= M 






Is 



M 

ii y^' y‘ 

^J=l 



Is 



Other posterior summau'ies caui be computed similarly from the replicated 
saimples, i and H having been selected to achieve "convergence'' to 
"smooth" estimates. 

4.4 A NUMERICAL EXAMPLE 

To illustrate the Gibbs sampler technique, we shall reauialyze the data 
set given by Kaplan and Meier (1958). The data consist of deaths 

occurring at .8, 3.1, 5.4 and 9.2 months and losses occxirring at 1.0, 
2.7, 7.0 and 12.1 months. For comparison purposes, we consider the 
SEime prior sjjecif ications used by Susarla and Van Ryzln (1976) in their 



with 4> 



Bayesian reanalysis of the data. That is, “ 1 “ 

.12 and N = 4,8, and 16. 

To apply the Gibbs saonpler approach, we divide the positive real line 
into the following intervals: (0, .8 ], (.8 , .8], (.8, 1], (1, 2.7], 

(2.7, 3.l”], (3.r, 3.1], (3.1, 5.4~], (5.4”, 5.4], (5.4, 7], 

(7, 9.2”], (9.2”, 9.2], (9.2, 12.1], and (12.1, «). We label these 

intervals by (0, t^], (t^, t^] ^^i2' *13^’ ®1* ®2’ ‘ ‘ 

respectively, denote the probabilities assigned to the Intervals. 



The likelihood of 0 is prop>ortional to 

U9) - . Sg 6j3) x 



(6g 9,3)(6 io 9,3>9j3 ■ 

-9t -ft 

Let ttj = W(e - e ), so that the prior distribution of 0 is 



n(0) = C 



13 a-1 

n 



j=i 



where C is the normalizing constant. 

Note that 0_, 0_, 0„, 0,, euid 0,_ in the likelihood combine simply with 

^ b o 11 IJ 

the corresponding 0 variables in the prior distribution, so that the 



parameters 0_, 0 



0 , 



2 ' 6 ’ 8 ’ 



0 , 



and 0^2 are each up-dated by 1 in the 

posterior distribution. Therefore, we need only introduce three Z 

variables for the incomplete data, namely, Z = (Z. ,, Z_ Z,_ ,), 

1 41 Ol IJ, 1 

We 



^2 ■ ^^ 52 ’ ^62 ^ 13 , 2 ^' ^3 



^^10,3* ^11,3' ^12,3’ ^13,3^' 



then sample Zj, for J - 1, 2, and 3, from the appropriate multinomial 
distribution with sample size 1 and rescaled probabilities. 



To estimate the survival function at t 



f 



we accumulate the 0j for i > 



J. For 



t between tj and interpolation formula 



that connects 



the survival function at the two end points according to the prior 
shap>e can be used. Tables 1 euid 2 exhibit the Gibbs sampler results 
for the survival function evaluated at tj with M - 1000 and H = 4000, 

both with i = 10. The exact Bayes solutions given by Susarla and Van 
Ryzin are also listed for compairison. The tables show that the Gibbs 
sampler results for H *= 1000 are already very accurate in approxima- 



ting the exact Bayes rules. SimilaLT results hold for N = 16. For 
further illustration of the Gibbs seunpler methodology, see Kuo (1991), 
who reanalyses data from Turnbull (1974). 



Table 1: Gibbs Approximation to the Bayes Estimates for Af = 4 



Statistics \ age(t) 


.8 


.8 


1 


2.7 


3. 1 


3. 1 


with H = 1000 


.970 


.886 


.879 


.819 


.805 


.702 


with M = 4000 


.970 


.886 


.879 


.819 


.805 


.701 


Exact Bayes 


.970 


.886 


.879 


.819 


.805 


.701 


Statistics \ age(t) 


5.4” 


5.4 


7 


9.2" 


9.2 


12. 1 



P^ with M = 1000 


.632 


.529 


.491 


.437 


.305 


.253 


P^ with H = 4000 ' . 


.632 


.529 


.491 


.438 


.307 


.256 


Exact Bayes 


.632 


.528 


.490 


.438 


.306 


.255 



Table 2: Gibbs Approximation 


to the Bayes 


Estimates 


for N 


= 8 


Statistics \ age(t) 


.8- 


.8 


1 


2.7 


3. r 


3. 1 


Pj with M = 1000 


.954 


.892 


.881 


.792 


.773 


.698 


P^ with H » 4000 


.954 


.892 


.881 


.792 


.773 


.700 


Exact Bayes 


.954 


.892 


.881 


.793 


.773 


.699 


Statistics \ age(t) 


5.4” 


5.4 


7 


9.2” 


9.2 


12. 1 


Pj with W = 1000 


.600 


.527 


.474 


.405 


.316 


.249 


Pj with M = 4000 


.602 


.529 


.474 


.405 


.318 


.250 


Exact Bayes 


.602 


.528 


.474 


.405 


.318 


.250 



ACKNOWLEDGEMENTS 



The work of the first author is supported by NSF grant DMS9008021 and 
the Naval Postgraduate School, Monterey, CA 93943, U. S. A. Both 
authors have benefited from discussions with A. E. Gelfauid. 



REFERENCES 

ANTONIAK, C. E. , (1974). "Mixtures of Dirichlet processes with applica- 
tions to Bayesian nonparametric problems". Arm. Statist., 2, 

1152-74. 

DELLAPORTAS, P. and SMITH, A. F. M. , (1991). "Bayesian inference for 

Generalised Linezir and Proportional Heizards Models via Gibbs 
Saunpling". Applied Statistics. To appeau'. 

FARD, N. S. and DIETRICH, D. L. , (1987). "A Bayesiaui note on 

reliability growth during a development testing prograjmne". 

I. E. E. E. Trans. Rel. , R-36, 568-572. 

FERGUSON, T. S. , (1973). "A Bayesian analysis of some non- 

parajnetric problems". Ann. Statist., 1, 209-30. 

f 

FERGUSON, T. S. and PHADIA, E. G. , (1979). "Bayesian nonparajnetric 

estimation bsised on censored data". Ann. Statist., 7, 163-86. 

GELFAND, A. E. and SMITH, A. F. M. , (1990). "Sampling-baised approaches 
to calculating marginal densities". J. Am. Statist. Assoc., 85, 
398-409. 

GELFAND, A. E. , HILLS, S. E. , RACINE-POON, A. and SMITH, A. F. M. , 
(1990). Illustration of Bayesiam inference in normal, data models 
using Gibbs saunpling. J. Amer. Statist. Assoc., 85, 972-985. 

GELFAND, A. E. and SMITH. A. F. M. (1991). "Gibbs sampling for 

marginal posterior expectations". Commun. in Statist. To app>ear. 

GELFAND, A. E. , SMITH, A. F. M. and LEE, T-M. , (1991). "Bayesian 
amalysis of constrained parameter and truncated data problems using 

Gibbs saunpling". J. Am. Statist. Assoc. To appear. 

GEMAN, S. and GEMAN, D. , (1984). "Stochastic relaxation, Gibbs 

distributions and the Bayesiaui restoration of images". I. E. E. E. 
Trans, on Pattern Analysis and Machine Intelligence, 6, 721-741. 

GILKS, W. R. and WILD, P. ,(1991). "Adaptive rejection saunpling for 
Gibbs Sampling". Applied Statistics. To appear. 

KAPLAN, E. L. and MEIER, P. , (1958). "Nonpau-^unetric estimation from 

incomplete observations." J. Am. Statist. Assoc., 53, 457-81. 



KUO, L. (1991). "Sampling based approach to computing nonparametric 
Bayes iain estimators with doubly censored data". Computing Science 
and Statistics: Proceedings of the 23rd Symposium on the Interface, 

SCHMEE, J. EUid HAHN, G. J. , (1979). "A simple method for regression 

analysis with censored data". Technometrics, 21, 417-432. 

SMITH, A. F. M. , (1977). "A Bayesiem note on reliability growth during 
a development testing programme". I. £. £. £. Trans. Rel. , R-26, 

346-347. 

SMITH, A. F. M. , SKENE, A. M. , SHAW, J. E. H. and NAYLOR, J. C. , 
(1987). "Progress with numerical and Graphical Methods for Bayes ism 
Statistics". The Statistician, 38, 75-82. 

SUSARLA, V. and VAN RYZIN, J. , (1976). "Nonpsuametric estimation of 

survival cxirves from incomplete observations. " J. Am. Statist. 
Assoc., 71, 897-902. 

TANNER, M. and WONG, W. , (1987). "The calculation of posterior 

distributions by data augmentation", (with discussion). J. Am. 

Statist. Assoc., 81, 82-86. 

TURNBULL, B. W. , (1974). "Honpsmametric estimation of a survivorship 

function with doubly censored data. " J. Am Statist. Assoc. , 69, 

169-73. 



(FILE: ASBCSK.CHI (JUMC /ppiJ) 



INITIAL DISTRIBUTION LIST 



1. Library (Code 0142) 2 

Naval Postgraduate School 

Monterey, CA 93943-5000 

2. Defense Technical Information Center 2 

Cameron Station 

Alexandria, VA 22314 

3. Office of Research Administration (Code 012) 1 

Naval Postgraduate School 

Monterey, CA 93943-5000 

4. Prof. Peter Purdue 1 

Code OR-Pd 

Naval Postgraduate School 
Monterey, CA 93943-5000 

5. Department of Operations Research (Code 55) 1 

Naval Postgraduate School 

Monterey, CA 93943-5000 

6. Prof. James D. Esary 1 

Code OR/Ey 

Naval Postgraduate School 
Monterey, CA 93943-5000 

7. Prof. Donald Gaver, Code OR-Gv 1 

Naval Postgraduate School 

Monterey, CA 93943-5000 

8. Prof. Patricia Jacobs 1 

Code OR/Jc 

Naval Postgraduate School 
Monterey, CA 93943-5000 

9. Prof. Harold J. Larson 1 

Code OR/La 

Naval Postgraduate School 
Monterey, CA 93943-5000 



10. Prof. P. A. W. Lewis 1 

Code OR/Lw 

Naval Postgraduate School 
Monterey, CA 93943-5000 

11. Prof. Paul R. Milch 1 

Code OR/Mh 

Naval Postgraduate School 
Monterey, CA 93943-5000 

12. Prof. Robert R. Read 1 

Code OR/Re 

Naval Postgraduate School 
Monterey, CA 93943-5000 

13. Prof. Lyn R. Whitaker 1 

Code OR/Wh 

Naval Postgraduate School 
Monterey, CA 93943-5000 

14. Prof. W. Max Woods 1 

Code OR/Wo 

Naval Postgraduate School 
Monterey, CA 93943-5000 

15. Prof. Thomas Ferguson 1 

Math Department 

University of California 
Los Angeles, CA 90024 

16. Prof. Robert Jennrich 1 

Math Department 

University of California 
Los Angeles, CA 90024 

17. Prof. Keh-Chau Li 1 

Math Department 

University of California 
Los Angeles, CA 90024 

18. Prof. Donald Ylvisker 1 

Math Department 

University of California 
Los Angeles, CA 90024 



19. Prof. Barry Arnold 1 

Statistics Department 

University of California 
Riverside, CA 92521 

20. Prof. Keh-Shin Lii 1 

Statistics Department 

University of California 
Riverside, CA 92521 

21. Prof. James Press 1 

Statistics Department 

University of California 
Riverside, CA 92521 

22. Prof. David Strauss 1 

Statistics Department 

University of California 
Riverside, CA 92521 

23. Prof. Joseph R. Gani 1 

Mathematics Department 

University of California 
Santa Barbara, CA 93106 

24. Prof. S. R. Jammalamadaka 1 

Statistics Department 

University of California 
Santa Barbara, CA 931067 

25. Prof. John Hsu 1 

Statistics Department 

University of California 
Santa Barbara, CA 931067 

26. Prof. Milton Sobel 1 

Statistics Department 

University of California 
Santa Barbara, CA 931067 



1 



27. Prof. Bradley Efron... 

Statistics Dept. 

Sequoia Hall 
Stanford University 
Stanford, CA 94305 

28. Prof. Iain Johnstone 1 

Statistics Department 

Stanford University 
Stanford, CA 94305 

29. Prof. Herbert Solomon 1 

Statistics Department 

Stanford University 
Stanford, CA 94305 

30. Prof. Byron Brown 1 

Biostatistics Department 

Stanford University 
Stanford, CA 94305 



31. Prof. George Roussas 1 

Statistics Division 

University of California 
Davis, CA 95616 

32. Prof. P. K. Bhattacharya 1 

Statistics Division 

University of California 
Davis, CA 95616 

33. Prof. Frank Samaniego 1 



Statistics Division 
University of California 
Davis, CA 95616 

34. Prof. Jane-Ling Wang 1 

Statistics Division 
University of California 
Davis, CA 95616 



1 



35. Dr. Jack C. Lee 

Bell Communications Research 
Morris Research and Engineering Center 
435 South Street, 2Q-374 
Morristown, NJ 07960 



36. Prof. Malay Ghosh 1 

Statistics Department 

University of Florida 
Cainvesville, FL 32611 

37. Prof. Glen Meeden 1 



Theoretical Statistics Department 
University of Minnesota 
270 Vincent 
206 Church St. SE 
Minneapolis, MN 55455 



38. Prof. James Zidek 1 

Statistics Department 

University of British Columbia 
Vancouver, V6T 1W5 CANADA 

39. Prof. Bruce Turnbull 1 



Operations Research and Industrial Engineering 
Cornell University 
227 FfTC 

Ithaca, NY 14853-3801 



40. Center for Naval Analyses 1 

4401 Ford Avenue 

Alexandria, VA 22302-0268 

41. Dr. Terry Speed 1 

Statistics Department 

University of California 
Berkeley, CA 94720 

42. Dr. P. Welch 1 

IBM Research Laboratory 



Yorktown Heights, NY 10598 



1 



43. Dr. Neil Gerr 

Office of Naval Research 
Arlington, VA 22217 



44. Dr. J. Abrahams, Code 1111, Room 607 1 

Mathematical Sciences Division, Office of Naval Research 

800 North Quincy Street 
Arlington, VA 22217-5000 

45. Prof. H. Chemoff 1 



Department of Statistics 
Harvard University 
1 Oxford Street 
Cambridge, MA 02138 

46. Prof. Carl N. Morris 1 

Statistics Department 

Harvard University 
1 Oxford St. 

Cambridge, MA 02138 

47. Prof. Tom A. Louis 1 

School of Public Health 

University of Minnesota 
Mayo Bldg. A460 
Minneapolis, MN 55455 



48. Prof. William Jewell 1 

Operations Research Department 

University of California, Berkeley 
Berkeley, CA 94720 

49. Dr. S. R. Dalai 1 

Bellcore 

445 South Street 
Morristown, NJ 07962-1910 

50. Dr. Michael P. Cohen 1 



National Center for Education Statistics 
555 New Jersey Ave., NW 
Washington, DC 20208-5654 



1 



51. Prof. Lynn Kuo 

Code OR-Ku 

Naval Postgraduate School 
Monterey, CA 93943-5000 



3 2768 00347504 7 



